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Abstract 


We report on a study of color coherence effects in pp collisions at a center 
of mass energy = 1.8 TeV. The data were collected with the D0 detector 
during the 1992-1993 run of the Fermilab Tevatron Collider. We observe the 
presence of initial-to-final state color interference with the spatial correlations 
between soft and hard jets in multijet events in the central and in forward 
pseudorapidity regions. The results are compared to Monte Carlo simulations 
with different color coherence implementations and to the predictions of 0(a 3 ) 
QCD calculations. 
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Color coherence phenomena in the final state have been very well established in e + e 
annihilations H, in what has been termed the “string” 0 or “drag” [§ effect. Particle 
production in the region between the quark and antiquark jets in e + e _ —> qqg events is 
suppressed. In perturbative quantum chromodynamics (QCD) such effects arise from inter¬ 
ference between the soft gluons radiated from the q, q, and g. While quantum mechanical 
interference effects are expected in QCD, it is important to investigate whether such effects 
survive the nonperturbative hadronization process, for a variety of reactions over a broad 
kinematic range, as predicted by local parton-hadron duality [§. 

The study of coherence effects in hadron-hadron collisions is considerably more subtle 
than that in e + e~ annihilations due to the presence of colored constituents in both the ini¬ 
tial and final states. The structure of multijet events in hard processes is influenced by the 
underlying color configurations at short distances. During a hard interaction, color is trans¬ 
ferred from one parton to another and the color-connected partons act as color antennae. 
Examples of color flow diagrams are shown in Fig. 1 for qq and qg scattering. In Fig. la 
(qq) the color system in which interference occurs is entirely between the initial and final 
states, whereas in Fig. lb (qg) interference also occurs both in the initial and in final states 
due to their explicit color connection. Gluon radiation from the incoming and outgoing 
partons forms jets of hadrons around the direction of these colored emitters. The soft gluon 
radiation pattern accompanying any hard partonic system can be represented, to leading 
order in 1 /N c where N c is the number of colors, as a sum of contributions corresponding to 
the color-connected partons. Within the perturbative calculations, this is a direct conse¬ 
quence of interferences between the radiation of various color emitters, resulting in the QCD 
coherence effects IBHTO 


Color coherence, which results in a suppression of soft gluon radiation in the partonic 
cascade in certain regions of phase space, can be approximated by Angular Ordering (AO). 
For the case of outgoing partons, AO requires that the emission angles of soft gluons decrease 
monotonically as the partonic cascade evolves away from the hard process. The radiation is 
confined to a cone centered on the direction of one parton, and is bounded by the direction 
of its color-connected partner. Outside this region the interference of different emission 
diagrams becomes destructive and the azimuthally integrated amplitude vanishes to leading 
order. Conversely the emission angles increase for the incoming partons as the process 
develops from the initial hadrons to the hard subprocess. Monte Carlo simulations including 
coherence effects probabilistically by means of AO are available for both initial and final state 
evolutions. Conventional perturbative QCD calculations taken to sufficiently high order 
should in principle incorporate these effects to any accuracy. Use of the latter approach is 
limited, however, due to the current lack of such calculations beyond order a^. 

Evidence for color coherence effects between initial and final states in pp interactions has 
been previously published [jlTJ for the central pseudorapidity region. In this paper we report 
on a new study of initial-to-final state color coherence phenomena in pp interactions which 
also extends the measurements into the untested forward region. We compare our central 
and forward measurements to the predictions of O(a^) QCD matrix element calculations 
as well as to several parton-shower Monte Carlo (MC) simulations. This is the first direct 
comparison of such measurements to analytic 0(a^) QCD predictions. 

0 


The D0 detector is described in detail elsewhere 


This analysis uses the uranium 
liquid-argon sampling calorimeter to measure jet final states. The D0 calorimeter has fine 
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segmentation in both azimuthal angle 0, and in pseudorapidity rj = — ln[tan(0/2)], where 6 
is the polar angle of the jet with respect to the proton beam. It has hermetic coverage for 
\rj 1 <4 with fractional transverse energy Et resolution of rv./ 80%/ \J '-Er(GeV) for jets. 

The data sample for this analysis [|T3[], representing an integrated luminosity of 8 pb ~ 1 , 
was collected during the 1992-1993 Tevatron Collider run. Events were selected using a three 
level trigger. The first level required a coincidence of two scintillator hodoscopes located on 
either side of the interaction region, to ensure an inelastic collision. The next stage required 
the transverse energy of at least three calorimeter towers (0.2 x 0.2 in A rj x A0) to exceed 
a 7 GeV threshold in the region \r)\ < 3.2. The surviving events were analyzed by an online 
processor farm where jets were reconstructed using a simplified version of the final jet Ending 
algorithm and an event was recorded to tape if it had a jet with Et > 85 GeV. 

Jets were reconstructed offline using an iterative fixed-cone clustering algorithm with 
cone radius R = \J(Arj) 2 + (A0) 2 = 0.5. Spurious jets from isolated noisy calorimeter cells 
and accelerator losses were eliminated by loose cuts on the jet shape. The E T of each jet was 
corrected for offsets due to the underlying event, multiple pp interactions, and noise; out- 
of-cone showering; and detector energy response as determined from the missing transverse 
energy balance of photon-jets events [14|. Events were required to have a measured vertex 
within 50 cm of the detector center. 

The remaining events were required to have three or more reconstructed jets. Jets were 
ordered in Et and labeled Et\ > Et 2 > Et 3 . We required Et\ > 115 GeV to avoid any 
bias introduced by the trigger threshold, and the third jet to have Et > 15 GeV. Color 
coherence effects are measured with the angular distribution in (rj, 0) space of the softer 
third jet around the second jet. The polar variables R and 0 = tan -1 ( Slgn ^ 3 ^ 32 ) were 
used to locate the third jet in a search disk of 0.6 < R < | around the second jet (Fig. 2). 
Here, ip and 0; are the pseudorapidity and azimuthal angle of the i th jet, A rj 32 = r/ 3 — rj 2 , 
and A032 = 03 — 02- We define 0 = 0 to point towards the beam nearest to the second 
jet, and 0 = n to point towards the farther beam. We study the interference effects in 
regions \r} 2 \ < 0.7 and 0.7 < l^l < 1.5. The pseudorapidity of the leading jet was not 
constrained, but the first and second jets were required to be in opposite 0 hemispheres, i.e. 
| < A02i < . After all selection criteria, a sample of 9,048 (5,776) events remains in the 

central (forward) region. 

The measured angular distributions (Fig. 3) are compared to the predictions of sev¬ 
eral MG simulations that differ in their implementation of color coherence. We employ the 
parton-shower MC event generators ISAJET 7.13 [p~5[ , HERWIG 5.8 |T7)j|, and PYTHIA 5.7 
a and a partonic event generator, JETRAD 1.2 [|18| to calculate the O(al) QCD predic¬ 
tions. The ISAJET generator uses an independent shower development model without any 
color coherence effects. Both HERWIG and PYTHIA incorporate initial and final state color 
interference effects by means of the AO approximation of the parton cascades. In PYTHIA, 
the AO constraint can be turned off. HERWIG and PYTHIA each employ a phenomenological 
model to describe the hadronization process. HERWIG uses the cluster hadronization model 
and PYTHIA implements the Lund string fragmentation model, which are both supported 
by the observations of color coherence phenomena in e + e“ annihilations. 

The shower-based MC simulations were performed at the particle (hadron) level after 
the nonperturbative hadronization process, whereas the JETRAD predictions were at the 
parton level. Detector rj and energy resolution effects were included in all predictions. The 
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generated events were subsequently processed using the same criteria employed for analyzing 
the data. 

The [3 distributions, normalized to the total number of events, for the data and HERWIG 
and ISAJET predictions are shown in Fig. 3 for both central and forward regions. In each 
case the data peak near (3 = tt and this enhancement is more dramatic for higher 77 . The 
shape of these distributions is sensitive not only to the process dynamics but also to phase 
space effects resulting from our jet and event selection criteria. For both regions, HERWIG 
is in good agreement with the data, whereas ISAJET shows systematic deviations from the 
observed distributions. 

The ratios of the observed f3 distributions relative to the MC predictions for both r 7 
regions are shown in Fig. 4. The data show a clear enhancement of events compared to 
ISAJET near the event plane (i.e., the plane defined by the directions of the second jet and 
the beam axis, f3 — 0, tt) and a depletion in the transverse plane ((3 — |). This is consistent 
with the expectation from initial-to-hnal state color interference that the rate of soft jet 
emission around the event plane be enhanced with respect to the transverse plane. The 
PYTHIA predictions include string fragmentation. Without AO the PYTHIA distributions 
are significantly different from the data, while with AO turned on there is much better 
agreement, although there are still some residual differences in the “near beam” region. 
In addition, from the H ^^ IG and JE trad @ distributions, we conclude that both the AO 
approximation and 0(al) QCD describe the coherence effects seen in data in both r] regions. 

The main sources of uncertainty on the data (3 distributions are summarized in Table |. 
Since we report event normalized distributions, any possible uncertainty on quantities that 
affect the overall rate of events is minimized. For both // regions the statistical and systematic 
uncertainties are comparable. Sources of systematic uncertainty arise from the jet energy 
calibration, a possible rj dependence of the jet reconstruction efficiency, and small // biases 
caused by the jet reconstruction algorithm. 

We have also examined three main sources of systematic uncertainty related to the MC 
predictions. Varying the jet energy resolution within its measured uncertainties resulted 
in changes to the MC (3 distributions of less than 2.6% (2.8%) in the central (forward) 
region. The effect of using different parton distribution functions (CTEQ2ML, CTEQ2MF, 
and CTEQ2MS [|Hj) was examined with JETRAD and found to produce variations of less 
than 2.0% in the central and forward regions. The JETRAD f3 distributions varied by less 
than 1.3% (2.5%) in the central (forward) region when the renormalization and factorization 
scales varied from E? i /2 to 2E T \. 

Table |TJ shows the y 2 values of fits to the various Mo ^ Garlo ratios with a constant for 
the central and forward regions as well as for the combined sample with |?/ 2 | < 1-5- All 
uncorrelated systematic uncertainties (i.e., energy scale and 77 bias corrections) in the data 
were added in quadrature with the statistical uncertainties and were included in the cal¬ 
culation of y 2 . The uncorrelated systematic uncertainty due to jet energy resolution was 
also added in quadrature with the statistical uncertainties of the MC predictions. From 
the y 2 values we conclude that HERWIG and JETRAD agree best with our data for both 
7 ] regions, whereas ISAJET clearly disagrees with the data. I 11 addition, from the data to 
PYTHIA comparisons we conclude that for the process under study, string fragmentation 
alone cannot accommodate the effects seen in the data. The AO approximation is an ele¬ 
ment of parton-shower event generators that needs to be included if color coherence effects 
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are to be modeled successfully. 

The data (3 distributions were also compared to the predictions of HERWIG at the parton 
level (before the nonperturbative hadronization stage). From the comparison of the H ^^ IG 
ratios at the parton and particle level we conclude that the hadronization effects as modeled 
by HERWIG are negligible and do not influence our results. 

We processed a limited HERWIG event sample through a full GEANT based detector 
simulation J2(J to investigate any possible detector effects not accounted for in the jet energy 
and r) resolutions. From the comparison of the /3 distributions at the particle and detector 
levels such residual detector effects were found to be negligible. 

We have presented a study of color coherent radiation in multijet events in pp collisions. 
We have measured the spatial correlations between the second and third leading Et jets in 
the central and in forward pseudorapidity regions. Comparisons of the data distributions 
with various color coherence implementations demonstrate a strong presence of initial-to- 
final state interference. Parton shower MC simulations that implement color interference 
by means of angular ordering reproduce the data angular distributions well, with HERWIG 
giving the best representation. Striking differences are found between the data and MC 
models that do not incorporate color coherence effects in the parton shower development. 
Our results also indicate that coherence effects as predicted by a 2 —> 3 partonic level 
calculation provide a good representation of our data, giving additional evidence supporting 
the validity of the local parton-hadron duality hypothesis. 

We express our deep appreciation to V. Khoze and T. Sjostrand for numerous valu¬ 
able discussions on this work. We also thank T. Sjostrand for helping us with the color 
coherence implementations in PYTHIA. We thank the staffs at Fermilab and collaborating 
institutions for their contributions to this work, and acknowledge support from the De¬ 
partment of Energy and National Science Foundation (U.S.A.), Commissariat a L’Energie 
Atomique (France), State Committee for Science and Technology and Ministry for Atomic 
Energy (Russia), CNPq (Brazil), Departments of Atomic Energy and Science and Education 
(India), Colciencias (Colombia), CONACyT (Mexico), Ministry of Education and KOSEF 
(Korea), and CONICET and UBACyT (Argentina). 
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TABLES 


TABLE I. 

Major uncertainties in the data (3 distributions. 


Source of uncertainty 

M < 0.7 

0.7 < \rj 2 \ < 1.5 

Jet energy scale corrections 

2 .2% 

3.4% 

Jet 7 ] bias correction 

1 .2% 

1.7% 

Jet reconstruction efficiency 



Near beam 

+1.5% 

+5.0% 

Far beam 

-0.5% 

-2.0% 

Statistical error 

3.2% 

4.3% 


TABLE II. x 2 values for 8 degrees of freedom of fits to the various Mon teCario ra L° s °f P distribu¬ 
tions with a constant for the central, forward, and combined samples. Statistical and uncorrelated 
systematic uncertainties were included. 

Event 


X 2 


Generator 

7?2 < 0.7 

0.7 < M < 1.5 

M < 1-5 

ISAJET 

48.0 

24.9 

74.4 

HERWIG 

5.1 

7.7 

7.2 

JETRAD 

7.9 

7.9 

7.2 

PYTHIA 




String fragmentation, but no AO 

77.2 

62.9 

126.0 

String fragmentation, and AO 

12.7 

11.3 

19.3 
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FIGURE CAPTIONS 


Figure 1 Color flow diagrams for (a) qq and (b) qg scattering. 


Figure 2 Three jet event topology illustrating the search disk (shaded area) for the softer 
third jet around the second jet. 


Figure 3 Comparisons of the data /3 distributions to the predictions of ISAJET and HERWIG 
for (a), (b) central region and (c), (d) forward region. The error bars include statistical errors 
only. 


Figure 4 Ratio of /3 distributions between data and the predictions of: (a) ISAJET, (b) 
PYTHIA with AO off, (c) PYTHIA with AO on, (d) HERWIG, (e) JETRAD for the central 
region; and (f)-(j) for the forward region respectively. The error bars include statistical and 
uncorrelated systematic uncertainties. 
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Figure 1 
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Figure 2 
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Figure 4 






